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ABSTRACT 

Methods  and  algorithms  for  navigation  systems  with  body-mounted 
inertial  instruments  have  been  studied  by  various  research  groups  throughout 
the  country.  The  high  iteration  rates  for  coordinate  transformation  in  a 
P3ERSHING  II  missile  and  the  complexity  of  the  navigation  solutions  make  a 
general  purpose  computer  or  special  purpose  computer  unfeasible  for  this 
application.  A  scheme  to  use  a  small  Digital  Differential  Analyzer  as  the 
front  end  of  a  general  purpose  navigation  computer  to  perform  the  coordinate 
transformation  is  feasible. 

This  report  describes  a  unique  mechanization  of  a  Digital  Differential 
Analyzer  to  transform  the  outputs  from  inertial  instruments  from  a  body-fixed 
coordinate  system  to  a  stabilized  coordinate  system.  The  design  is  specially 
adaptable  to  various  levels  of  logic  density  using  large-scale  integrated 
circuits  because  it  has  few  different  logic  circuits  used  repeatedly  and  performs 
serial  computation  that  requires  minimum  input  and  output  leads.  Unique 
features  of  the  Digital  Differential  Analyzer  are  that  it  solves  a  3  x  3  direction 
cosine  matrix  by.  merely  circulating  data  around  a  pair  of  registers  and  adders 
for  each  matrix  element  and  uses  a  two-bit  adder  to  add,  subtract,  and 

multiply  by  ( 2) n. 
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1.  Introduction 

Various  research  groups  throughout  the  country  have  made 
feasibility  studies  and  developed  algorithms  for  the  computations  required  for 
a  system  with  body-mounted  inertial  sensors  (strapdown).  The  Systems 
Branch  of  the  Army  Inertial  Guidance  and  Control  Laboratory  and  Center  has 
also  made  such  studies  and  established  requirements,  computational  speeds 
and  accuracies,  and  selected  algorithms  for  a  computer  to  perform  those 
computations. 

The  high  iteration  rates  for  computational  accuracy  of  those  algorithms 
in  a  PERSHING  II  application  plus  the  complexity  of  computation  of  the  naviga¬ 
tion  equations  make  the  use  of  a  general  purpose  computer  or  special  purpose 
computer  unfeasible.  A  scheme  to  use  a  small,  fast,  Digital  Differential 
Analyzer  which  employs  fast,  low-powered,  compact,  large-scale  integrated 
circuits  to  compute  the  direction  cosine  algorithms  was  selected.  This 
Digital  Differential  Analyzer  will  transform  the  outputs  of  inertial  sensors 
from  a  body-fixed  coordinate  system  to  a  stabilized  coordinate  system  and 
feed  them  to  a  simple  navigation  computer  for  solution  of  the  balance  of  the 
guidance  and  control  problem. 

This  report  describes  the  mechanization  of  a  Digital  Differential 
Analyzer  designed  to  be  built  from  large-scale  integrated  circuits  to  perform 
the  coordinate  ti  ansformation.  More  precisely,  this  Digital  Differential 
Analyzer  was  designed  to  transform  the  outputs  of  body-mounted  gyroscopes 
and  accelerometers  for  missiles  up  to  a  PERSHING  II  class  in  a  strapdown 
application  into  equivalent  outputs  from  a  stabilized  gimbal  system.  Beyond 
thaf  navigation  and  guidance  problem  remains  unchanged. 

.he  algorithms  and  organization  for  this  transformation  are  described 
in  another  report. 1  The  scheme  to  mechanize  those  algorithms  is  also 
described  in  another  report. 2  This  latter  report  describes  the  various  schemes 
considered  and  tradeoffs  effected,  the  type  of  outputs  from  the  inertial  instru¬ 
ments,  and  the  speeds  required  to  maintain  computational  accuracy  required 
for  such  a  missile  system. 


1U.  S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama,  Digital 
Differential  Analyzer  Organization  for  an  Analytic  Platform  Transform  Com" 
puter  by  Gerald  Scheiman,  January  1967,  Report  No.,  RG-TR-67-2] 

2U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama,  Pershing 
Semi-Annual  Report,  Guidance  Systems  Branch  by  M.  C.  Jones,  et  al. ,  8  Febru¬ 
ary  1968,  Report  No.  ITN-RGG-18. 
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The  flow  and  sequence  of  data  representing  the  two  updatings  of  the 
algorithms  given  below  are  shown  in  Figure  1.  The  complete  computer  has 
three  rows  operating  simultaneously.  A  more  detailed  description  of  each 
updating  appears  later  in  this  report.  The  appendix  gives  the  flow  of  data. 


TIMING  PERIOD  0 


TIMING  PERIOD  2 


AY  ACCUMULATOR 
FOR  ROW  ONE 


O  ADDER 
CZ3  REGISTER 


UPDATING  1 

NOTE:  ROWS  TWO  AND  THREE  ARE  IDENTICAL 


FIGURE  1.  FLOW  OF  DATA  FOR  TWO  UPDATINGS  OF  DIRECTION 
COSINES  FOR  ROW  ONE 


TIMING  PERIOD  5 
UPDATING  2 


NOTE:  SEE  APPENDIX  FOR  FLOW  OF  DATA  IN  COMPLETE 
COMPUTER  FOR  TWO  UPDATINGS 

FIGURE  1.  FLOW  OF  DATA  FOR  TWO  UPDATINGS  OF  DIRECTION 
COSINES  FOR  ROW  ONE  (Concluded) 


The  algorithms 
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where 


i  =  1,  2,  and  3  (for  the  three  strings) 

C.  (K)  =  value  of  cosine  matrix  element 

1  j  J. 

th 

(K)  --  the  K  cycling  of  the  total  algorithm,  updatings  1  and  2 
A0V,  A0V,  A 0  =  incremental  outputs  from  the  body-mounted  gyros. 

Algorithms  are  for  solving  3x3  direction  cosine  matrix.  Only  one  of  the 
three  strings  is  shown  in  two  successive  updatings.  The  other  two  strings 
are  identical  and  the  updating  sequence  repeats. 


2.  Description 

The  block  diagram  (Figure  2)  represents  the  complete  Digital 
Differential  Analyzer  and  interface  with  the  test  computer  and  the  navigation 
computer.  There  are  nine  matrix  elements.  The  hardware  for  each  of  the 
matrix  elements  is  identical,  each  being  made  up  of  one  pair  of  adders  and 
one  pair  of  16-bit  registers.  They  are  arranged  in  three  identical  columns 
and  three  rows.  At  any  given  time,  the  operation  of  every  row  is  identical 
and  the  columns  differ  only  by  the  number  of  clock  pulses  or  by  whether  they 
are  outputing  or  inputing.  The  sequencer  and  1/ O  provide  the  control  functions. 
The  following  discussions  will  be  limited  to  only  one  row  since  the  other  two 
are  physically  identical.  The  clocking  for  any  column  is  the  same  for  any 
timing  period  and  is  applied  to  each  pair  of  adders  simultaneously. 

The  computer  uses  a  2?s  complement  number  system  in  which  negative 
numbers  (sign  bit  =  1)  are  the  2's  complement  of  positive  binary  numbers 
(sign  bit  =  0) .  This  makes  it  possible  to  use  a  simple  adder  to  perform 
addition  and  subtraction  and  automatically  get  the  correct  sign  in  the  result. 

In  subtraction  the  computer  automatically  forms  the  2*s  complement  of  the 
subtrahend  and  adds.  A  short  explanation  of  this  system  is  given  in  Table  I. 


3.  Sequencing  and  Operations 
a.  Organization 

The  studies  and  trade-offs  in  the  previously  mentioned  report3 
established  the  word  length  for  direction  cosine  to  be  16  bits  plus  sign  and  the 

3 Jones,  loc.  cit. 
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FIGURE  2.  BLOCK  DIAGRAM  OF  THE  DIGITAL  DIFFERENTIAL 
ANALYZER  SHOWING  PRESET  VALUES  OF  DIRECTIONAL 
COSINES  IN  PRFSET  POSITION 


TABLE  I.  EQUIVALENTS  OF  BINARY  AND  DECIMAL  SYSTEMS 


— 

BINARY 

DECIMAL 

BINARY 

DECIMAL 

01111 

15 

10001 

-15 

OHIO 

14 

10010 

-14 

01101 

13 

10011 

-13 

01100 

12 

10100 

-12 

01011 

11 

10101 

-11 

01010 

10 

10110 

-10 

01001 

9 

10111 

-  9 

01000 

S 

liOOO 

-  8 

00111 

7 

11001 

-  7 

00110 

6 

11010 

-  6 

00101 

5 

11011 

-  5 

00100 

4 

11100 

-  4 

00011 

3 

11101 

-  3 

00010 

2 

11110 

-  2 

00001 

1 

11111 

-  1 

00000 

0 

NOTE:  To  form  2’s  complement  (negative  form)  of  a  binary  number, 
start  from  the  least  significant  end;  copy  all  0  bits  and  the  first 
bit;  then  complement  all  the  rest. 

Example: 

Binary  10  01010 

2’s  complement  10110 

An  example  of  2’s  complement  arithmetic  may  be  shown: 

a.  (+5)  +  (+7)  =  +12 

+  +5  00101 

+7  00111 

+12  01100 

b.  (  +  5)  +  (-7)  =  -2 

+5  =  00101 
+  -7  =  11001 
-2  =  11110 


(+5)  - 

(-7)  =  +5 

+7  =  +12 

<->  -+75 

00101 

+  5  = 

00101 

11001 

+  +7  = 

00111 

+12  =  01100 


(-5) 

— 

(+7)  =  -5 

-7  =  -12 

(-) 

-6 

11011 

-5  =  11011 

+7 

00111 

-7  =  11001 

-12  =  10100 

basic  clock  to  be  1.  2  mHz.  The  direction  cosines  of  the  A 0  column  partially 

A 

increment  the  A0V  and  A 0  columns  in  one  word  time;  the  A 0V  column  par- 

tially  increments  the  A 0  and  A 0  columns  and  the  A0  column  partially  incre- 

x  z  z 

ments  the  A0V  and  A0  columns  (Algorithms  list,  paragraph  1) .  Complete 
X  Y 

updating  is  obtained  in  three  word  times,  after  which  the  sequence  is  reversed 
to  increase  numerical  accuracy  of  the  computation  when  limit  cycle  angular 
inputs  are  received.4  The  direction  cosines  from  each  column  are  combined 
with  appropriate  corresponding  AV^  (body  axis)  to  convert  them  to  AVg  (space 
axis) . 


b.  Scaling  and  Updating  of  Direction  Cosines 


The  gyros  are  scaled  so  that  an  angular  increment  of  A0  will 
output  one  pulse  and  increment  the  direction  cosine  by  a  maximum  of  one  bit 
in  the  bit  16  position.  Since  the  direction  cosines  update  each  other,  the 
matrix  elements  were  so  sized  that  each  direction  cosine  word  is  actually  31 
bits  plus  sign  in  length,  and  these  are  held  in  pairs  of  16 -bit  registers  that  are 
clocked  simultaneously  during  operation.  During  the  first  word  period  the 
most  significant  part  of  the  A0_,  column  is  added  to  the  least  significant  part  of 

the  A 0V  and  A0„  columns.  The  overflow  resulting  from  such  addition  is  the 

X  Z 

carry  left  in  the  respective  adder  and  is  automatically  added  to  the  most 
4Scheiman,  loc.  cit. 
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significant  part  of  the  A0  and  A 9  columns  to  increment  them  partially.  This 

Y  Z 

is  a  very  convenient  arrangement,  since  the  adder  is  the  only  place  where  a 

bit  can  be  temporarily  stored.  The  next  timing  period  of  the  register  will  shift 

and  pick  up  the  stored  carry.  Similarly,  during  the  second  word  time,  the 

A0  ma  "ix  element  will  partially  update  the  A0  and  A0  matrix  elements,  and 
Y  “  X  Z 

during  the  third  word  time,  the  A 0  matrix  element  matrix  will  partially  update 

z 

the  A0  and  A 0  matrices  element,  thus  completing  one  updating  of  direction 

-A  X 

cosines  in  one  row.  The  same  occurs  simultaneously  in  the  other  two  rows. 

The  maximum  accelerometer  output  occurs  less  than  half  as  often  as  the  gyros 
and  are  sampled  to  be  combined  with  direction  cosines  every  other  updating 
period. 


c.  Sequencing 

The  updating  arrangement  is  straightforward,  but  mechanization 
is  complex  because  of  the  unequal  shifting  that  is  necessary  at  every  column 
during  every  wofd  time  and  the  change  in  scale  factor  of  the  gyros  for  outputs 
during  maximum  body  rates.  Figure  3  illustrates  how  the  most  significant 
pare  of  the  outputing  matrix  element  combines  with  the  least  significant  part 
of  the  matrix  element  which  is  partially  incremented.  The  sign  bit  is  added  to 
all  bits  until  it  coincides  with  the  sign  bit  of  the  matrix  element  being  updated. 
There  is  a  slight  overlap,  and  the  matrix  element  that  is  to  output  next  must 
stop  when  bit  16  is  the  next  bit  through  the  adder.  This  uneven  clocking  requires 
overlapping  to  complete  any  partial  updating.  The  overlapping  occurs  during 
the  next  word  time  when  the  next  partial  updating  takes  place.  Therefore,  it 
is  necessary  to  keep  track  of  A0,  sign  bit,  bit  31,  carry,  and  carry  reset.  The 
direction  cosine  at  the  top  is  updated  by  the  direction  cosine  at  the  bottom. 

The  pulse  rebalanced  gyros  being  considered  for  this  application  will 
double  the  scale  factor  when  body  rate  exceeds  a  certain  level  in  that  axis. 


MOST  SIGNIFICANT  PART  LEAST  SIGNIFICANT  PART 


31 

17 

16 

1 

S . ss 

S 

31 

16 

MOST  SIGNIFICANT  PART 


FIGURE  3.  COMPARISON  OF  THE  MOST  AND  LEAST  SIGNIFICANT 
PARTS  OF  THE  OUTPUTING  MATRIX  ELEMENT 
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Pulses  from  the  gyro  will  then  represent  2A0,  and  the  direction  cosines  must 
be  updated  twice  as  much  or  twice  as  fast  as  those  for  the  A 0  case.  In  most 
cases,  this  change  in  scale  factor  will  not  occur  in  all  three  u>.<»s  at  the  same 
time,  so  the  choice  of  doubling  iteration  rates  at  random  for  one,  two,  or  all 
three  axes  is  not  feasible.  This  is  analogous  to  shifting  gears  on  a  tricycle 
but  affecting  only  one,  two,  or  sometimes  all  three  wheels;  of  course,  the 
machine  could  be  made  to  operate  always  twice  as  fast  as  in  the  worst  case 
and  provide  logic  to  sample  gyro  pulses  every  other  time  for  the  A0  case. 

This  approach,  however,  would  unnecessarily  increase  logic  speed.  Logic 
for  multiplying,  even  by  two,  is  much  more  complex  than  the  plain  two-bit  adders 
being  used.  The  mechanization  implemented  for  the  2A0  was  to  double  the 
size  of  the  updating  direction  cosine  by  running  it  through  a  flip-flop,  causing 
a  one-clock  period  delay  as  it  entered  the  input  adder  of  the  matrix  element 
being  updated.  In  binary  numbers  this  is  the  same  as  shifting  the  whole  word 
one  bit  to  the  left,  which  is  alternately  the  same  as  multiplying  by  two.  The 
solution  was  simple,  but  logic  had  to  be  added  for  the  computer  to  recognize 
the  2A0  case  and  delay  the  word  without  changing  the  number  of  clock  pulses 
for  the  A0  case  and  retain  the  information  during  the  overlap  portion  of 
updating.  The  word  time  for  register  will  vary  from  15  to  17  clock 
pulses.  Two  successive  updatings  are  accomplished  in  six  timing  periods  (0  to 
5)  of  twenty  clock  pulses  (tO  to  tl9)  each.  The  sequencing  stops  each  register 
at  the  appropriate  time.  The  following  sequence  describes  what  takes  place 
in  each  matrix  element  of  each  row  during  every  timing  period.  It  is  based  on 
the  assumption  that  A0  pulse  is  present  in  the  matrix  unit  being  undated  as  will 
be  the  case  with  binary  pulse  rebalance  gyros.  The  effect  of  zero  angular  input 
to  the  gyro  will  be  that  the  gyro  outputs  and  algorithms  will  alternate  plus  and 
minus  for  a  net  direction  cosine  increment  of  zero. 


(1)  Timing  Period  £  (Figure  4) 


(a)  Overlap.  During  previous  timing  period  5  adder  A2 
outputs  bits  16  to  31  into  matrix  element  4  (ME4)  and  ME7  through  A3  and  A5, 
respectively.  The  sign  bit  was  not  added,  so  it  is  the  first  one  out  of  register 
(R2)  during  this  timing  period.  Therefore,  the  first  clock  period  will  allow 
the  sign  bit  to  go  into  A3  and  A5  to  set  the  sign  to  be  added  to  bits  17  plus  sign 
and  complete  the  partial  updating  started  during  5. 

(b)  Partial  Updating  This  Period.  Matrix  element  1  will 
receive  17  clock  pulses  and  output  bits  16  to  31  plus  sign  through  A1  to  partially 
increment  ME4  and  ME7.  Matrix  element  4  will  be  the  next  to  output  through 
A4,  so  it  will  receive  only  15  clock  pulses  and  thus  add  only  bits  16  to  30  from 
ME1  to  bits  1  to  15.  Matrix  element  7  receives  16  clock  pulses  because  it  will 


be  updated  through  A5  during  the  next  timing  period  and  must  have  bit  1  as  the 
next  bit  into  A5.  Bits  16  to  31  from  ME1  are  added  to  bits  1  to  16. 

During  this  riming  period,  bits  16  to  31  plus  sign  are  simultaneously 

outputed  to  the  AV  accumulator  where  they  are  combined  with  AV  to  form 

XB 

the  AV  component  of  AV  ready  to  be  sampled  by  the  navigation  computer. 

X  Xb 

(2)  Timing  Period  jL  (Figure  5) 

(a)  Overlap.  Bits  31  and  sign  from  ME1  into  ME4  and 
sign  bit  from  ME1  into  ME7  must  be  added  to  complete  the  partial  updating  of 
A0  and  A0  columns  from  timing  period  0.  Therefore,  the  first  two  clock 

X  £t 

pulses  of  this  period  will  add  bit  31  and  set  sign  in  A4  of  ME4  and  the  first 
clock  pulse  of  this  period  will  set  sign  in  A6  of  ME7. 

(b)  Partial  Updating  This  Period.  Matrix  element  1 
receives  18  clock  pulses  to  place  bit  1  as  the  next  one  into  A2  since  it  will  also 
receive  data  during  the  next  timing  period  2.  Bits  18  to  31  from  ME4  are 
added  to  bits  1  to  16  in  Al.  Matrix  element  4  receives  17  clock  pulses  to  out¬ 
put  through  A4  to  partially  increment  ME1  through  Al  and  ME7  through  A5. 
Matrix  element  7  receives  15  clock  pulses  to  stop  bit  16  just  before  A5  because 
it  will  increment  ME4  and  ME1  during  the  next  timing  cycle.  Bits  16  to  30 
from  A4  are  added  to  bit3  1  to  15  through'  A5. 


(3)  Timing  Period  2  (Figure  6) 

(a)  Overlap.  This  overlap  is  the  same  as  that  of  timing 
period  1  except  that  it  includes  bits  31  and  sign  from  ME4  into  ME7  and  sign 
from  ME  4  into  ME1. 

(b)  Partial  Updating  This  Period.  Matrix  element  1 
receives  16  clock  pulses  to  place  bit  one  as  the  next  bit  into  Al  since  it  will 
also  receive  data  during  the  next  timing  period  !3.  Bits  16  to  31  from  ME7  are 
added  to  bits  1  to  16  in  A4  of  ME4,  and  in  A2  of  ME1. 

During  this  timing  period,  bits  16  to  31  plus  sign  are  simultaneously 

outputed  to  the  AV  accumulator  where  they  are  combined  with  AV  to  form 
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the  AV  component  of  AV  . 
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(4)  Timing  Period  3  (Figure  7) 


(a)  Overlap,  The  first  clock  pulse  adds  the  sign  bit 
from  ME7  into  ME4  and  ME1  to  complete  the  partial  updating  of  A 0^  and  A 6^ 
from  timing  period  2. 


(b)  Partial  Updating  This  Period.  Matrix  element  1 
receives  16  clock  pulses  to  place  bit  1  as  the  next  bit  into  A2  since  it  will  also 
receive  data  during  the  next  timing  period  3.  Bits  16  to  31  from  ME7  are 
added  to  bits  1  to  16  of  ME1  through  Al.  Matrix  element  4  receives  15  clock 
pulses  to  place  bit  16  as  the  next  bit  into  A3.  Bits  16  to  30  from  ME7  are 
added  to  bits  1  to  15  of  ME4  through  A3.  Matrix  element  7  receives  17  clock 
pulses  to  output  bits  16  co  31  plus  sign  and  place  bit  1  as  the  next  bit  into  A6 
since  it  will  receive  data  during  the  next  timing  period. 

( 5)  Timing  Period  4  ( Figure  8) 


(a)  Overlap.  The  overlap  of  timing  period  4  is  the 
same  as  that  of  timing  period  1  except  that  it  is  bit  31  and  sign  from  ME7  into 
ME 4  and  sign  bit  from  ME7  into  ME1. 

(b)  Partial  Updating  This  Timing  Period.  Matrix 
element  1  receives  15  clock  pulses  to  place  bit  16  as  the  next  bit  into  A2  since 
it  will  output  data  during  the  next  timing  period.  Bits  16  to  30  from  ME4  are 
added  to  bits  1  to  15  of  ME1  through  A2.  Matrix  element  4  receives  17  clock 
pulses  to  output  bits  16  to  31  plus  sign  and  put  bit  1  as  the  next  bit  into  A3 
since  it  will  receive  data  during  the  next  timing  period.  Matrix  element  7 
receives  16  clock  pulses  to  make  bit  1  the  next  bit  into  A5  since  it  will  receive 
data  during  the  next  timing  period.  Bits  36  to  31  from  ME4  are  added  to  bits 

1  to  16  of  ME7  through  A6. 

During  this  timing  period,  bits  16  to  31  plus  sign  are  simultaneously 

outputed  from  ME4  to  the  AV  accumulator  where -they  are  combined  with  AVv_, 
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to  form  the  AV  component  of  AVvc,. 


( 6)  Timing  Period  5  ( Figure  9) 


(a)  Overlap.  The  first  clock  pulse  adds  the  sign  bit 

from  ME4  into  ME7  and  bit  31  and  sign  from  ME4  into  ME1  to  complete  the 

partial  updating  of  A 0  and  A 6  from  timing  period  4. 

X  z 


receives  16  clock  pulses  to  output  bits 


ib)  Updating  This  Timing  Period.  Matrix  element  1 


16  to  31  and  place  bit  16  as  the  next  bit 
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into  A1  since  it  will  also  output  during  the  next  timing  period.  Matrix  element 
4  receives  16  clock  pulses  to  place  bit  1  as  the  next  bit  into  A4  since  it  will 
also  receive  data  during  the  next  timing  period.  Bits  16  to  31  from  M£1  are 
added  to  bits  1  to  16  in  ME4  through  A3.  Matrix  element  7  receives  16  clock 
pulses  to  place  bit  1  as  the  next  bit  into  A6  since  it  will  also  receive  data 
during  the  next  timing  period.  Bits  16  to  31  from  ME1  are  added  to  bits  1  to 
16  in  ME7  through  A5. 

The  overlap  from  this  timing  period  was  described  above  with  timing 
period  (L  That  will  conclude  two  successive  updatings.  Operations  will  con¬ 
tinue  by  recycling  through  the  six  timing  periods. 


4.  Interface  With  Navigation  Computer 

The  direction  cosines  and  increments  of  velocity  generated  in  this 
Digital  Differential  Analyser  are  held  in  registers  during  operation.  Data  are 
sampled  by  the  navigation  computer  for  solving  navigation  equations  and  storage 
in  memory.  Logic  in  the  Digital  Differential  Analyzer  is  arranged  fr>r  auto¬ 
matic  presetting  and  readout  of  registers  by  the  navigation  computer. 

The  D84  PERSHING  ground  computer  will  be  used  for  laboratory  tests. 
However ,  any  comparable  general  purpose  computer  can  be  used.  The  Digital 
Differential  Analyzer  contains  a  buffer  register  for  parallel  input  or  output. 

For  presetting  or  register  readout  during  testing,  the  navigation  computer  or 
test  computer  need  only  be  programmed  to  select  a  matrix  element  line  for 
presetting  or  register  readout.  The  Digital  Differential  Analyzer  will  auto¬ 
matically  accept  incremental  gyro  and  accelerometer  outputs  and  output 
direction  cosines  and  increments  of  velocity  as  required  by  the  navigation 
computer. 


5.  Hardware  Design 
a.  Registers 

There  are  eighteen  serial  registers,  each  16  bits  long, and  six 
serial  17-bit  registers.  All  are  identical  with  output  parallel  of  the  first  two 
most  significant  bits  as  required.  The  eighteen  are  used  in  pairs  in  each  of 
the  nine  matrix  elements  and  the  six  are  used  in  pairs  in  each  of  the  AV 
accumulators.  All  registers  will  be  built  alike  and  the  sequences  and  external 
lines  used  will  adapt  them  to  the  matrix  location  or  AV  accumulator  where 
u-sed.  The  I/O  has  one  17-bit  register  with  capability  for  parallel  and  serial 
input  and  output. 
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b.  Adders 


There  are  eighteen  identical  two-bit  adders  used  in  pairs 
with  each  pair  of  the  above  registers  in  each  of  the  nine  matrix  elements.  The 
sequencer  and  input  leads  adapt  the  adder  design  for  use  in  any  matrix  location. 
With  the  2's  complement  system  and  minor  modifications,  this  adder  adds, 
subtracts,  and  multiplies  input  data  by  2.  The  AV  accumulator  also  uses  a 
two-bit  adder,  but  its  inputs  are  different  enough  to  justify  a  slightly  different 
design. 

c.  Synchronization 

The  basic  computer  timing  is  provided  by  a  1.  2-mHz 
clock,  a  recycling  counter  that  counts  clock  pulses  0  to  19  (20  clock  pulses) , 
and  a  recycling  counter  that  counts  20^  clock  pulses  and  goes  from  0  to  5 
(six  timing  periods) .  All  timing  functions  are  derived  from  combinations 
of  these  basic  counters  driven  by  the  clock. 

d.  Sequencer 

The  sequencer  is  shown  as  a  block  in  Figure  2;  however, 
it  is  made  up  of  several  circuits  driven  by  the  counter  or  clock  and  working  in 
various  combinations  to  sequence  each  matrix  element  through  updating  and 
overlap  cycles.  The  sequencing  of  the  three  matrix  elements  in  any  column  is 
identical  for  any  given  timing  period;  therefore,  one  set  of  circuits  makes  up 
the  sequencer  as  follows: 

( 1)  Sign  Set  and  Reset  and  Carry  Reset.  Updating  of  a 
direction  cosine  occurs  only  if  A0  occurred  and  each  A0  has  a  sign  to  indicate 
which  way  the  direction  cosine  should  be  updated.  The  algorithms  alternate 
sign  of  A 0  to  equalize  direction  cosine  updating  when  there  are  no  body  rates; 
therefore,  the  sign  of  A0  is  inverted  so  that  the  updating  of  the  direction  cosine 
is  accomplished  with  the  proper  sign  (+)  (-)  =  -,  (-)  (-)  ■=  +. 

As  pointed  out  before,  the  complete  direction  cosine  updating  (six  word 
times)  consists  of  two  three-word-time  updatings.  Therefore,  each  adder 
inputs  twice  during  the  complete  updating  and  always  at  the  same  timing  period 
of  the  six  timing  period  cycle.  During  the  input  cycle,  logic  circuitry  sets  or 
resets  a  flip-flop  in  the  adder  when  the  sign  bit  is  processed  through  the  adder. 
The  flip-flop  holds  the  sign  bit  through  the  overlap  until  the  sign  bit  of  the 
direction  cosine  is  processed.  The  same  logic  that  resets  the  sign  at  that 
point  also  resets  the  carry  flip-flop  in  the  adder. 
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(2)  Logic  for  Initialization  and  Sign  for  A 9  and  2A 9.  The 
sign  of  A0  is  held  in  a  series  of  six  flip-flops  that  are  updated  once  every  six 
timing  periods.  There  is  a  seventh  flip-flop  in  this  logic  with  a  load  or  lift-off 
state.  The  load  position  inhibits  all  the  adders  while  the  Digital  Differential 
Analyzer  is  being  preset.  The  lift-off  position  sets  one  flip-flop  for  each  of  the 
three  axes  and  resets  one  flip-flop  for  each  of  the  three  axes.  With  no  gyro 
motion,  the  binary  pulse  rebalance  gvros  will  alternate  ±A 9  and  produce  a  zero 
net  updating  of  direction  cosines.  Gyro  motions  will  cause  predominance  in 
plus  or  minus,  causing  updating  in  that  direction.  The  logic  for  A 9  is  six 
flip-flops  initialized  in  the  Ad  position.  The  gyro  inputs  are  sampled  every 
six  timing  periods  and  will  change  to  2A 9  when  body  rates  exceed  specified 
limits. 


6.  Input/Output 

The  I/O  consists  of  nine  flip-flops  to  select  a  matrix  element  for 
presetting,  nine  flip-flops  to  select  a  matrix  element  for  readout,  and  three 
flip-flops  to  select  the  AV  register  for  readout  of  the  accumulated  velocity. 

A  buffer  register  holds  the  data  that  are  being  preset  or  read  out  of  the  Digital 
Differential  Analyzer.  The  buffer  register  is  enabled  only  when  one  of  the  21 
flip-flops  has  been  selected.  Data  are  clocked  in  or  out  of  the  buffer  once 
with  exactly  17  clock  pulses  until  another  flip-flop  is  selected. 


7.  Hardware  Construction 

The  Digital  Differential  Analyzer  was  designed  to  be  built  with 
Targe-scale  integrated  (LSI)  circuits.  The  fact  that  all  the  matrix  elements  are 
identical,  except  for  slight  variations  in  input  signals,  offers  several  levels 
of  logic  density  concentration  that  can  be  implemented  with  LSI  circuits.  The 
design  has  been  partitioned  for  a  breadboard  of  one  row  plus  controls,  to  be 
built  out  of  T2L  integrated  circuits  components.  This  will  allow  laboratory 
testing  and  debugging  before  the  whole  design  Is  committed  to  any  particular 
LSI  circuit  approach.  While  the  breadboard  is  built  and  tested,  progress  of 
the  various  LSI  circuit  approaches  is  being  closely  monitored  to  effect  a 
better  selection.  The  LSI  circuit  approach  to  be  selected  will  not  necessarily 
be  the  one  with  the  highest  gate  concentration.  It  is  mere  probable  that  it  will 
be  one  with  a  reasonably  high  gate  concentration,  inherent  high  reliability, 
low  cost,  and  realistic  approach  to  logistics. 
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8.  Nuclear  Effect  Considerations 

The  shift  registers  making  a  high  portion  of  the  computer  are, 
unfortunately,  the  most  sensitive  components  to  nuclear  radiation.  The  gamma 
and  X-ray  radiation  will  very  likely  cause  change  of  state.  Techniques  for 
hardening  are  available  but  not  simple.  Since  the  navigation  computer  will  be 
using  all  the  current  data  from  the  Digital  Differential  Analyzer,  it  will  be 
more  feasible  to  use  the  memory  in  the  navigation  computer  for  storage  and 
take  advantige  of  the  hardening  that  must  be  provided  there.  After  the  critical 
period  of  a  nuclear  blast,  the  information  can  be  reset  correctly  in  the  regis¬ 
ters.  This  is  the  most  logical  approach  since  the  navigation  computer  is  to 
provide  whatever  permanent  storage  the  Digital  Differential  Analyzer  requires 
and  since  it  also  reduces  the  problem  of  nuclear  hardening  to  a  single  solution. 


9.  Conclusion 

The  iteration  rates  to  solve  direction  cosines  for  a  PERSHING  II 
class  missile  in  a  strapdown  application  and  the  solution  of  the  navigation  prob¬ 
lem  would  be  prohibitive  for  a  state  of  the  art  general  purpose  machine. 

However,  the  high  speed,  low  power,  and  compact  features  of  LSI  circuits 
offer  an  excellent  approach  to  build  a  small  serial  Digital  Differential  Analyzer 
to  solve  the  direction  cosine  matrix  to  serve  as  inputs  to  the  navigation  computer. 
The  design  of  such  a  Digital  Differential  Analyzer  described  in  this  report  lends 
itself  idenlly  to  various  levels  of  LSI  circuit  construction  because  of  its  few 
different  logic  circuits  used  repeatedly  as  common  building  blocks,  and  its 
serial  operation  that  requires  minimum  of  input  and  output  leads. 

This  design  is  unique  in  that  the  3x3  direction  cosine  matrix  is  con¬ 
tinually  updated  by  merely  circulating  each  word  around  a  pair  of  adders  and 
registers  comprising  each  matrix  element  and  simultaneously  accepting  or 
outputing  data  to  matrix  elements  in  the  same  string  coincident  with  sampling 
of  gyro  inputs.  The  use  of  a  two-bit  adder  to  add,  subtract,  and  multiply 
throughout  the  Digital  Differential  Analyzer  is  another  advantage  of  the  design. 
The  capability  to  multiply,  by  merely  adding  a  flip-flop  at  the  input  to  the  adder, 
to  accommodate  random  changes  in  scale  factor  of  g>ro  inputs,  is  an  extremely 
simple  mechanization  when  compared  to  the  computer  complexity  that  would  be 
required  by  the  alternate  solutions  cortsidered.  The  alternate  solution  would 
have  made  the  mechanization  of  this  scheme  extremely  difficult  and  maybe 
prohibitive  with  present  day  hardware. 

The  mechanization  of  the  scheme,  as  unique  as  it  may  appear,  did  not 
get  by  without  complications.  It  would  have  been  very  desirable  to  complete 
every  incremental  updating  within  a  single  timing  period,  without  overlap.  The 
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need  for  overlapping  to  complete  the  partial  incrementation  became  apparent 
early  in  the  design  by  simulating  the  flow  of  data  on  paper.  Overlapping  to 
complete  partial  updatings  created  complications  because  it  is  not  always  the 
same.  It  requires  the  computer  to  store  the  sign  and  A 0  or  2A0  of  the  current 
and  previous  complete  updating  periods  (three  timing  periods) ;  it  also  requires 
proper  sequencing  for  each  of  the  columns.  Fortunately,  the  logic  to  accom¬ 
plish  all  the  sequencing  is  fairly  simple  and  is  a  one-time  application  in  the 
computer. 

Nuclear  hardening  is  recognized  as  a  problem  with  this  computer.  The 
approach  selected  is  to  keep  the  current  data  from  this  computer  stored  in  the 
navigation  computer  and  take  advantage  of  the  hardening  provided  for  the 
navigation  computer  rather  than  solving  the  same  problem  for  the  two  com¬ 
puters. 

The  complexities  if  the  computer  are  more  in  following  the  sequence  of 
operations  than  hardware  complexity.  The  design  can  be  implemented  using 
any  of  the  LSI  circuit  approaches  being  developed.  It  is  also  feasible  to  build 
each  row  with  a  different  LSI  circuit  approach  to  more  conclusively  evaluate 
the  various  more  promising  LSI  circuit  approaches. 
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Appendix 

FLOW  OF  DATA  IN  COMPLETE  COMPUTER  FUNCTION 
FOR  TWO  UPDATINGS 
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This  report  describes  a  unique  mephanization  of  a  Digital  Differential  Analyzer  to  trans¬ 
form  the  outputs  from  inertial  instruments  from  a  body-fixed  coordinate  system  to  a  stabilized 
coordinate  system.  The  design  is  specially  adaptable  to  various  levels  of  logic  density  using 
large-scale  integrated  circuits  because  it  has  few  different  logic  circuits  used  repeatedly  and 
performs  serial  computation  that  requires  minimum  input  and  output  leads.  Unique  features 
of  the  Digital  Differential  Analyzer  are  that  it  solves  a  3  x  3  direction  cosine  matrix  by  merely 
circulating  data  around  a  pair  of  registers  and  adders  for  each  matrix  element  and  uses  a  two- 
bit  adder  to  add,  subtract,  and  multiply  by  (2)n. 
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